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1. Introduction 

Amplitudes for hard scattering processes in Quantum Chromodynamics (QCD) are of basic 
importance both for theory and phenomenology and precision predictions for them must include 
higher-order quantum corrections. An important aspect in explicit computations are the singu- 
lar limits of amplitudes at higher loops. Here, one has to consider two types of limits, soft and 
collinear, related to the emission of gluons with vanishing energy and to collinear partem radiation 
off massless hard partons, respectively. 

For massless QCD amplitudes the corresponding singularities are regularized by working in d 
dimensions and appear as explicit poles in (d — 4). Typically two powers in l/(d— 4) are generated 
per loop. When massive particles are involved, some of the collinear singularities are screened by 
the parton masses, which gives rise to large logarithmically enhanced contributions. In both cases 
the structure of the singularities can be understood from the factorization property of QCD and 
it can be predicted to all orders based on a small number of perturbatively calculable anomalous 
dimensions. Moreover, factorization gives rise to an extremely simple universal multiplicative re- 
lation between a massless amplitude and its massive version in the limit when the parton masses 
are small with respect to all other kinematical invariants. This relation can be employed to derive 
virtual QCD corrections at higher loops including all logarithms in the heavy quark mass as well as 
all constant (mass-independent) contributions. We demonstrate the predictive power of factoriza- 
tion with recent results for the two-loop QCD amplitudes for heavy-quark production in hadronic 
collisions. 



2. Factorization of QCD amplitudes 

We are interested in general 2^n scattering processes of partons p\ 

p: -» P3~\ VPn+2- (2.1) 
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,v = -"p I {^l^m/l^c,}, -^2,a s (At 2 ),£ ) , (2.2) 



2]. In ^-dimensions we can write Eq. (2.2) as a product of functions ^ m °\ y^ m °^ and Jf?&\ 
>(m =o) = ^=0) ^ jflt(M 2 )jfi ^ ^{k^M^e) W, (2-3) 



The corresponding scattering amplitude ^ p depends on the set of fixed external momenta 
masses {m,} and color quantum numbers {q}, 

f 
A* 1 

as well as on the strong coupling constant a s , the renormalization scale ]X and the parameter e of 
dimensional regularization, d = 4 — 2e. Also, we denote explicitly the hard scale Q of the process 
typically related to the center-of-mass energy, e.g. Q = y 7 ? with s = (k\ +^) 2 - 

Let us briefly recall the factorization of on-shell amplitudes for massless partonic processes [1, 

r (m=0) ^{n 
■ p ) ^ p 

(,,.,-(>> . , 7 «;«-()! / 0" N / /1_ 

VF" / V H-~ J 

( m= n) 

where we use matrix notation suppressing the color indices. The jet function ^ depends 
only on the external partons. It collects all collinearly sensitive contributions and is color-diagonal. 
Coherent soft radiation arising from the overall color flow is summarized by the soft function 
J^p m ~°\ which is a matrix in color space. The short-distance dynamics of the hard scattering is 
described by the (infrared finite) hard function J4? p , which to leading order is simply proportional 
to the Born amplitude. 

The factorization formula (23) organizes the singularity structure of any massless QCD am- 
plitude. From the operator definitions for the functions ^ p and S^" 1-0 ^ and the corresponding 
renormalization group properties one derives evolution equations. The solution of the latter gives 
rise to an all-order exponentiation in terms of well-known anomalous dimensions. As an upshot, 
all 1/e terms related to the emission of gluons with vanishing energy and to collinear parton radi- 
ation off massless hard partons, respectively, exponentiate to all orders in perturbation theory, see 
e.g. [2,3]. 

Most important for our considerations is the jet function, which contains all collinear contri- 
butions from the external partons. It is therefore of the form 

a (w=o) = n 4r 0) = n Hr 0) )"> (2 - 4) 

z'G {all legs} ie {all legs} 

where i = q,g for quarks and gluons. ^h, is the individual jet function of each external parton, 



which by definition, one identifies with the (gauge invariant) form factor of a quark or a gluon, 
^ Of course, the latter function is well-known in QCD, see e.g. [4,5]. 

When masses are introduced the picture described above gets modified. However, the basic 
factorization of the QCD amplitude in jet, soft and hard function from Eq. (2.3) can be retained. 
With the exception of contributions related to heavy quark loops (see below) in the presence of a 
hard scale Q we can write for the partonic process (2.1 ) 

p) W = A m) f^WA(M 2 ),e)y p W (M,^,a s (At 2 ), £ ) », (2.5) 

where all non-trivial mass dependence enters in the functions and J?p , while power sup- 

pressed terms in the parton masses are neglected in Jfc. The jet function for massive partons can 
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be defined in complete analogy to Eq. ( |2.4| ), i.e. we identify w i tn tne massive form factor 



■ This guarantees exponentiation with largely the same anomalous dimensions as in the mass- 
less case. Also the soft anomalous dimensions which govern the soft function have a smooth limit 
for vanishing parton masses. 

i one can deduce a remarkably simple rela- 



In summary, from comparison of Eqs. ( [Q| ) and (£ 
tion between a massless and a massive amplitude in the small-mass limit. Thus, QCD factorization 
provides us with [6] 



,(») 



n 

(6 {all legs} 



7 (m\0) 
1*1 



2 X Jtv™ . 



where we have again suppressed the color indices, ^p'" -0 ^ and the corresponding massive am- 
plitude in the small mass limit m 2 <C Q 2 are multiplicatively related by a universal process 
independent function Z^°\ With the definitions for the jet functions in Eq. (p~4|), it is directly 
given in terms of the respective form factors, 



(2.6) 



,(m|0) 



m 



•>2 2 



<2 2 



-l 



(2.7) 



where the index i denotes the (massive) parton and a s = a s (/i 2 ). The process-independence is 



manifest in Eq. ( |2.7[ ), because Z^"' ^ is only a function of the (process-independent) ratio of scales 
jl 2 /m 2 . The (process-dependent) scale Q cancels completely between the massive and the massless 
form factors. 



Eq. (2.6) can be used to predict any massive amplitude from the known massless one, which 



is a great advantage in practice, as the latter is much easier to compute. Moreover, Eq. (2.6) in- 
cludes not only the singular terms and the logarithms in the massive amplitude but extends even 
to the mass-independent constant contributions. However, there is one important side condition on 
Eqs. (^^) and ( |2.7[ ), concerning all terms proportional to the number of heavy quarks n h . These 
two-loop contributions are excluded explicitly from the definition of Z^T . In order to incorpo- 
rate them additional process dependent terms appear as has been shown e.g. in QED for Bhabha 
scattering [7]. 



3. Hadro-production of heavy quarks 



As an application of the formalism developed, we consider the pair-production of heavy quarks 
in the ^-annihilation and the gluon fusion channel, 



q : 

g : 



q(ki)+q(k 2 ) 
g{ki)+g(k 2 ) 



Q(k 3 ,m) + Q(k 4 ,m) 
Q(k 3 ,m) + Q(k 4 ,m) 



(3.1) 



where kj denote the on-shell parton momenta and m the mass of the heavy quark, thus k\ 
and k 2 



k 2 = m 2 . Energy-momentum conservation implies k^ +k% = k 1 ^ +k% and we consider the 



scattering amplitude for the processes ( p.lD in QCD perturbation theory, 

2 



>(m) 



(0)\(w) 



+ 



(2)v(n 



+ 



a. 



(3.2) 



4 



Massive QCD amplitudes at higher orders 



S. Moch 



which defines the series expansion in the strong coupling a s = a s (/i 2 ) and ji. is the renormalization 
scale. As usual the MS-scheme for the coupling constant renormalization is employed and the mass 
m is taken to be the pole mass. 

It is convenient to define the function £/ p (e,m,s,t,jj,) for the squared amplitudes summed over 
spins and colors as 



£|^ p | 2 = s/ p (e,m,s,t,n), p = q,g- 

/ v is a function of the Mandelstam variables s, t and u given by 



(3.3) 



s= (P1+P2) 



t = (Pi ~Pi) 



m 



u = (P1-P4) 



m 



(3.4) 



and has a perturbative expansion similar to Eq. 



). 



2™2 



M m ) 



271 J' 



/6,(m) 



ofe Y 

271 J 



/p 8,(m) + 



(3.5) 



In terms of the amplitudes the expansion coefficients in Eq. (3.5) may be expressed as 



4- {m) 


= {JtfX 


<) 


(m) ^ 






= (Jtf\ 


4 l) ) 


(m, + (.^p (,) |- 


#p (0) ) H , 


4^ 




4 l) ) 


W + (^p (0) |. 


4 (2) }W + (4(V P (0) } (ffl 



(3.6) 
(3.7) 
(3.8) 



where the results for have been presented e.g. in [8, 9] and the so-called loop-by-loop 

contribution in i/ 8 '( m ) can be found in [10], both results with the complete dependence on the 
heavy-quark mass. The new contribution is the real part of {^p^\^p^)^ up to powers 0{m) in 
the heavy-quark mass [11, 12]. 

In order to obtain (^ p (0) \Jf^ 2) ) {n 

r (m\0) 



from Eq. ( |2.6| ), we have to construct the appropriate func- 



tions Zj.j from the on-shell heavy-quark form factor and the massless on-shell ones, all results 
being known [4,5, 12, 13] to sufficient orders in a s and powers of e. An explicit expression for 

2 



+ 2n M + \2n) 



a; 



(3.9) 



up to two loops is known [6]. As mentioned above, the definition ( p.9D accounts in particular for 
all fermionic terms except for those linear in n h . The leading ny terms ~ (n^a s )" for the process 
gg — > QQ in Eq. (3.1) can also be predicted, where we denote the total number of flavors with 
rif, which is the sum of n t light and n h heavy quarks. Keeping only terms quadratic in n h and/or 
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n h + n l one has up to two loops, 



where 



7 (m|0) 

A) - 



(l^/Osy (2) + 3) 



,(1) 

1«] 



with zP9 ~ n h known from [6, 12]. Note that is also equal to the <^(a s ) term in the gluon wave 
function renormalization constant Z3. 



2k M ^\27i) bl 



Hi, X M 



(3.10) 



(3.11) 
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Exploiting the predictive power of the relation Eq. (2.6) and applying it to the processes (3.1) 
we get 



2Re(X, (0) |A (2) ^ (m) 



+Z 



(i), 



A(m=0) 



2Re(X, (0) |A (2) ) (m=0) 



+ 2z£U 4,{m=0) + 



K 1 x«/°) + 



2Re(^ g (0) |^ g (2) )( m) 
,/ z (i) +z (i) 



2Re(.# g (0) |^ g (2) )( m=0) 

,6,(m=0) _ / (2) 7 (2) 7 (1) 7 (1) 
g + 1 \ L \Q] +%] + Z [g] Z [g] 



m 



/ 4.(m=0) 



+ 



+ 



(3.12) 



(3.13) 



m 



which assumes the hierarchy of scales m <C i.e. we neglect terms 6(m). Eqs. ( |3.12| ) and 



( J3- 1 3[ ) predict the complete real part of the squared amplitudes {jfcf' |^# p )W except (as indi- 



cated) for those terms, which are linear in n h . The real part of the two-loop massless amplitudes 
(^#p°^|^#p 2 ^ m=0 ^ are computed in [14, 15]. The finite remainders of the latter agree with the cor- 
responding terms constructed from the two-loop helicity amplitudes calculated of [16, 17] after the 
infrared subtraction procedure is performed. 

In order to arrive at a complete prediction for (^#p°^|^#p 2 ')( m ) including all heavy-quark loop 
corrections, the factorization approach in Eqs. ( J3.12 ), ( 3.13| ) is supplemented by a direct calcula- 
tion of all necessary massive Feynman diagrams as an expansion in the small mass. The advantage 



of this approach is an independent check of Eqs. (3.12) and (3.13) as well as of the correspond- 
ing massless results. It relies on the reduction of integrals to a set of masters with the Laporta 
algorithm [18], the subsequent construction of Mellin-Barnes representations for all the integrals, 
see e.g. [19-23] and the summation of series representations [24] or the application of the PSLQ 
algorithm [25]. 

We are now able to give the result for the interference of the two-loop and Born amplitude for 
the scattering processes ([0|). For a SU(./V)-gauge theory with N denoting the number of colors, 
one has C A = N, C F = (N 2 — \)/2N and T F = 1/2 and, as mentioned above, the total number of 
flavors tif = rij + n h is the sum of n f light and n h heavy quarks. Thus, for qq — > QQ and gg — > QQ 
we have 



2Re 



(0), sAT-Urn) 



2(N 2 - 1 ) f N 2 A q + B q + — C q + Nn,D qJ + Nn h D q>h 



+ ^E qJ + ^E qM + (n l + n h ) 2 F q 



2Re 



Co) i ^( 2 )\H 



(N 



1) lN\+NB g + ic g + ^D g +N\E g j +N 2 n h E g j l 



(3.14) 



(3.15) 



+n l F g j + n h F g ^ h + — G gJ + — j G gift 

n 2 

+Nn l 2 H gJ +Nn,n h H gMl +Nn 2 H gih + j^I g . 



n l n h n h 



and all explicit expressions can be found in [1 1, 12]. 

As explained above, the factorization approach provides results for all coefficients except the 



terms linear in n h . These are D q j, and E q j, in Eq. (3.14) and Egj,, F g j t and G g ,h in Eq. (3.15), which 
have been obtained from a direct calculation of the massive loop integrals as briefly sketched above. 



Furthermore, to have an independent cross check of the factorization formulae (3.12) and (3.13) 
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the direct computation of Feynman diagrams has also been extended to the coefficients A q , D q j, 



E q j and F q in Eq. ( |3.14| ) and to A g , E g j, H g j, H g ik, Hg,h, Ig,i, Ig,ih, and I g j r in Eq. ( p.l5| ). Of course, 
for the coefficients tested we have found full agreement between both methods. 



4. Conclusions 

We have presented a discussion of the singular behavior of on-shell QCD amplitudes with mas- 
sive particles at higher orders and we have emphasized the strong similarities between scattering 
amplitudes with massless and massive partons in the small mass limit. In this regime, factorization 
relates the two amplitudes multiplicatively by a universal process independent function Z^°\ The 
present results for amplitudes generalize the massless formulae known previously [1,2] and they 
extend one-loop massive results [26] to all orders. For cross sections the analogous property [27] 
is now known through next-to-next-to-leading order [28,29]. 

We have illustrated the predictive power of the factorization ansatz with new results for heavy- 
quark hadro-production at two-loops in QCD [11, 12]. The results include all logarithms in the 
mass m as well as the constant terms and they can be used as a strong check of any future complete 
calculation. Moreover, if combined with the threshold behavior of the amplitude, which is not 
known at present they can serve as a well founded basis for quantitative predictions for, say, top 
production at LHC to next-to-next-to-leading order. In order to obtain physical cross sections, 
the virtual amplitudes considered here have to be combined with the corresponding real emission 
contributions, of course. 

In a similar spirit, one can also derive higher order QCD corrections to new heavy (colored) 
particles, like squark or gluino production in supersymmetric extensions of the Standard Model. 
Also the approach taken here may prove useful in the future in formulating subtraction schemes 
with massive partons for the real emission contributions beyond one loop. 



References 



[1] S. Catani, Phys. Lett. B427 (1998) 161-171, [hep-ph/9802439] 



[2] G. Sterman and M. E. Tejeda-Yeomans, Phys. Lett. B552 (2003) 48-56, [hep-ph/0210130| 



[3] S. Mert Aybat, L. J. Dixon, and G. Sterman, Phys. Rev. D74 (2006) 074004, [hep-ph/0607309| 



[4] S. Moch, J. A. M. Vermaseren, and A. Vogt, JHEP 08 (2005) 049, [hep-ph/0507039|. 



[5] S. Moch, J. A. M. Vermaseren, and A. Vogt, Phys. Lett. B625 (2005) 245-252, [ nep-ph/0508055 1 



[6] A. Mitov andS. Moch, JHEP 05 (2007)001, [hep-ph/0 61214 9] 



[7] T. Becher and K. Melnikov, JHEP 06 (2007) 084, [arXiv : 0704 . 3582 [hep-ph] ]. 



[8] J. G. Korner and Z. Merebashvili, Phys. Rev. D66 (2002) 054023, [ |hep-ph/0207054| ]. 
[9] W. Bernreuther, A. Brandenburg, Z. G. Si, and P. Uwer, Nucl. Phys. B690 (2004) 81-137, 



[ Piep-ph/0403035Q . 



[10] J. G. Korner, Z. Merebashvili, and M. Rogal, Phys. Rev. D73 (2006) 034030, [ |hep-ph/0511264| | 

[11] M. Czakon, A. Mitov, and S. Moch, Phys. Lett. B651 (2007) 147-159, [ arXiv: 07 05 .1975 
[hep-ph] ]. 



7 



Massive QCD amplitudes at higher orders 



S. Moch 



[12] M. Czakon, A. Mitov, and S. Moch, [arXiv: 0707 . 4139 [hep-ph] ]. 



[13] W. Bernreuther et. at, Nucl. Phys. B706 (2005) 245-324, [ |hep-ph/0406046p . 

[14] C. Anastasiou, E. W. N. Glover, C. Oleari, and M. E. Tejeda-Yeomans, Nucl. Phys. B601 (2001) 
318-340, [ |hep-ph/0010212| ]. 

[15] C. Anastasiou, E. W. N. Glover, C. Oleari, and M. E. Tejeda-Yeomans, Nucl. Phys. B605 (2001) 
486-516, [|hep-ph/0101304|]. 



[16] Z. Bern, A. De Freitas, and L. J. Dixon, JHEP 06 (2003) 028, [hep-ph/0304168]. 



[17] A. De Freitas and Z. Bern, JHEP 09 (2004) 039, [hep-ph/0409007 [. 



[18] S. Laporta, Int. J. Mod. Phys. A15 (2000) 5087-5159, [nep-ph/0102033|. 



[19] V. A. Smirnov, Phys. Lett. B460 (1999) 397-404, [hep-ph/9905323]. 



[20] J. B. Tausk, Phys. Lett. B469 (1999) 225-234, [hep-ph/9909506). 



[21] M. Czakon, J. Gluza, and T. Riemann, Phys. Rev. D71 (2005) 073009, [hep-ph/0412164|. 



[22] M. Czakon, Comput. Phys. Commun. 175 (2006) 559-571, [hep-ph/0511200]. 



[23] M. Czakon, J. Gluza, and T. Riemann, Nucl. Phys. B751 (2006) 1-17, [hep-ph/0604101]. 



[24] S. Moch and P. Uwer, Comput. Phys. Commun. 174 (2006) 759-770, [ |math-ph/0508008| ]. 

[25] H. R. R Ferguson and D. H. Bailey, (see e.g. 

http : / /mathworld . wolfram . com/PSLQAlgorithm. html ) . 



[26] S. Catani, S. Dittmaier, and Z. Trocsanyi, Phys. Lett. B500 (2001) 149-160, [ |hep-ph/0011222| ]. 
[27] B. Mele and R Nason, Nucl. Phys. B361 (1991) 626-644. 



[28] K. Melnikov and A. Mitov, Phys. Rev. D70 (2004) 034027, [hep-ph/0404143]. 



[29] A. Mitov, Phys. Rev. D71 (2005) 054021, [hep-ph/0410205]. 



8 



